Bounds on the volume fraction of the 
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Abstract 

We study the size estimate problem for the two-phase shahow sheh 
equations in this paper. Our aim is to derive bounds on the volume 
fraction of each phase assuming that the material properties of the 
two phases are given. The approach in this paper is based on the 
translation method. One of the key steps is to connect the shallow 
shell equations to the thin plate equation. 



1 Introduction 

Let an elastic shallow shell consist of two materials. We can assume that one 
of them is an inclusion whose size is unknown. A fundamental question is to 
estimate the volume fraction of the inclusion by a noninvasive method such as 
measuring boundary responses. In this paper we want to prove some bounds 
on the volume fraction based on the translation method. To begin, we let 
f2 be a bounded simply connected domain in with smooth boundary dQ. 
For a shallow shell, its middle surface is described by {{xi,X2,e9{xi,X2)) : 
{xi, X2) G fl} for e > 0, where £ is a small parameter. Let u = {ui, U2, U3) = 
(u'jMa) : r2 — 7- represent the displacement vector of the middle surface. 
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Taking e — t- with appropriate scalings, Ciarlet and Miara [1] showed that 
u satisfies the following equations: 

fdiv = in il, 

[divdiv m-div(s^V^) = in 

where the constitutive laws are given as 

4/i 2 , 4A/i 2 w 

= 4^e^ + -i^(Tre^)/2, (1.2) 



A + 2/i 

^ (Vu + (Vu')^ + ® + (ye ® Vu^f) 



and X,fi are Lame coefficients. Here, a® b denotes the tensor product of 
vectors a and b and I2 is the 2x2 identity matrix. The Cauchy data on dfl 
associated with (1.1) is given by 

(u';n3,U3,„) (Dirichlet), 

(s^n; (div m — s^V6') ■ n + (mn ■ t)^t, mn ■ n) (Neumann), 

where n is the boundary normal vector field and t is the tangential vector 
field in the positive orientation. Also, we define g^n = ^9 ■ n and g^t = ^9 ■ t 
for any function g. Here and elsewhere in the paper Ca ■ b means (Ca) ■ b 
for any 2x2 matrix C. It is clear that (div m — s^V^) ■ n + (mn • t)^^ can 
be replaced by (div m) ■ n + (mn ■ t) ^ since s^n is already given. 
Assume now that the shell consists of two phases, i.e., 

A = AiXi + A2X2 and /x = niXi + Ai2Ai2, 

where 

Xj 



1 in phase j, 
otherwise. 



Here we are interested in the inverse problem of estimating the volume frac- 
tion of the first phase /i (or the second phase /2) by one pair of Cauchy data 
(1.3) when Ai, A2,/ii,/i2 are given. The approach taken here is based on the 
variational method, precisely, the translation method. On one hand, this pa- 
per extends results in [3] (for conductivity) and in [7] (for elasticity), which 
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used the translation method of Murat and Tartar [9, 10, 11] and Lurie and 
Cherkaev [5, 6], to the shallow shell equation. On the other hand, we would 
like to give another application of the authors' recent result on the equiva- 
lence of the inverse boundary problems for 2D elasticity and the thin plate 
[4]. Using the method in [4], the idea here is to transform the elasticity-like 
equation of (1-1) into a plate-like equation (4th order). At the same time, the 
Cauchy data for the new equation are completely determined by the original 
one. Therefore, the original problem is reduced to the size estimate problem 
for a decoupled system of 4th order equations, which is easier to deal with. 
For the shallow shell equations with an isotropic inhomogeneous medium in 
which the medium inside of the inclusion is not known, bounds on the volume 
fraction involving constants which are not explicitly given were derived in [2]. 
The method used in [2] is based on some quantitative uniqueness estimates 
for the shallow shell equations. 

This paper is organized as follows. In Section 2, we transform the first 
equation of (1.1) into a 4th order equation and discuss the corresponding 
Cauchy data for the new system. The inverse problem is then investigated 
in Section 3. 



2 Transformation of the equation and the Cauchy 
data correspondence 

It was known that the 2D elasticity equation and the thin plate equation are 
equivalent (see for example [8, Section 2.3]). The transformation of the first 
equation of (1.1) into a plate-like equation will follow the same idea. Denote 

^6 q6 
'11 ''12 

'12 ■'22 

The first equation of (1-1) is written as 

.0 , ^0 _ n 

'11,1 "T ■'12,2 — U, 
■^12,1 + -^22,2 ~ 0- 



Hence, there exists a function ip such that 

11 ■^12^ _ f '^,22 -1p,12 



'21 "^22 



-^,12 ^,11 



MV^^, (2.1) 
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where the fourth order tensor R is defined by 



= R^MRi 



for any 2x2 matrix M, where 

Ri 



r 
-1 oy 



The second equation of (1.2) can be written as 

= §s^ (2.2) 
where the fourth order tensor S is defined by 

§A = —A'^"" r(TrA)J2 (2.3) 

for any matrix A. On the other hand, from the form of (the third equation 
of (1.2)), we have 

_ 1 (V^ ® Vug + {V9 ® Vu-sf) = ^ (Vu' + (Vu')^) (2.4) 
and it follows that 

divdiv (Me^) - ^divdiv {R{Ve ® Vug) + M(V^ ® Vu^f) = 0. (2.5) 
Using (2.1) and (2.2) in (2.5) yields 

divdiv(LVV) - ^divdiv (M(V^ ® Vug) + M(V^ ® Vug)^) = 0, (2.6) 

where L = MSM. It is easily seen that L = §. Replacing the first equation of 
(1.1) by (2.6), (1.1) is transformed into 

f divdiv (LV^^;) - idivdiv (M(V^ ® Vua) + M(V^ ® Vug)^) = 0, 
|divdiv(MV2n3) - diviRV^ijjVe) = 0, ^ ' 



where 



We now discuss the Cauchy data corresponding to (2.7). The Cauchy 
data for (2.7) are given by 

div(LV^V^) ■ n - ^div (M(V^ ® Vug) + R(V^ ® Vug)^) ■ n 
+ (^(LVV)n-t- ^ (M(V^® Vms) +M(V^® Vms)^) n-t^ ; '''^'^^ 

(LVV) 

n n 

2 

and 

(div(MV\3) - MV^^V^) ■ n + ((MV^Ms)!! ■ t),^; (2.9) 
(MV^U3)n • n. 

The Cauchy data (2.9) correspond to the second equation of (2.7), which 
are aheady given in (1.3). Now we want to investigate (2.8). In view of 
authors' recent result [4], we can show that (2.8) are completely determined 
by (u', s^n) on dQ. In fact, by (2.1), one can determine ip and ip^n from s^n. 
On the other hand, let e = |(Vu' + (Vu')^), then we have 



Thus we obtain 
div RTeR 



U2,12 — \ui,22 — 1^2,12 \ (\U2,12 — \ui,22 



which implies 



^^^1,12 — 1^1,21 — 1^2,11 / \|lil,12 — 1^2,11^ 

({\ui,2-\u2,l),l\ 1 1 

R± Ll 1 , = R±V(-Mi,2 - -M2,l) 

V( 2^1,2 - 2^2,l),2/ ^ ^ 



(divRl^Rx) ■ n = (R^n) ■ V(^mi,2 - ^^2,1) 
= • V(-ni,2 - -U2,i) 

= (-M2,l - 2«l,2),i- 
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Combining (2.1), (2.2), and (2.4) yields that 

div(LV^^/') ■ n - ^div (M(V^ ® Vug) + M(V^ ® Vug)^) ■ n 

+ (^(LVV)n ■ t - ^ (M(V^ ® Vms) + M(V^ ® Vu^f) n ■ (2.10) 

= (^^2,1 - ^Mi,2),t + ((RleR^n) • t)^^ . 
Observe that |(m2,i — ""1,2) can be expressed as 

We also have 

(Rl^R^n) ■ t = f Rl ^"^'2 + '^""'A R^n 

From (2.10) we thus obtain 

div(LV^^^) ■ n - ^div (M(V^ ® Vug) + R(V^ ® Vug)^) ■ n 
+ |^(LV^V^)n • t - ^ (M(V^ ® Vna) + M(V^ ® Vug)^) n ■ t 

= -(t ■ Rl(Vu')t),i = -m±t) ■ (Vu')t),t = -(n ■ (Vu')t),i = -(n ■ u',),^. 

(2.11) 

We can also deduce that 

(LV^^)n n-- (R(V^ ® Vwg) + M(V^ ® VMa)^) n ■ n 
= n • (R^eRx)n = t ■ et = t ■ (Vu')t = t ■ u'^. 

Formulae (2.11) and (2.12) show that the Neumann data in (2.8) are com- 
pletely determined by u'^. 

It is helpful to take a closer look at Vug in the Neumann data of (2.8). It 
is obvious that VM3 (on dQ) is determined by U3 and u^^n- In other words, the 
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live ® vu^r ^ ^ ^''""^'^ "^'^"^^'^ 



values of (R(V^ ® VM3) + M( ® Vug)^) n-t and (M( ® Vwg) + R{V9 ® Vwg)^) n- 
n are given in terms of M3 and u^^n- It remains to simplify 

div (E(V^ ® Vms) + M(V^ ® Vu-sf) ■ n. 

Notice that 

R{V9 ® Vug) 

and 

Thus, we have that 

div(M(V^ ® Vug)) ■ n 

^ (6',2M3,2),l - (^,1%,2),2 \ ^ 

V-(^,2^3,l),l + (^,1^3,1)2; ' 

= (-^^21^3^ 12^3^1) ' ^ (^.2M3,21 - 0,lU3^22)ni + (-6',2M3,11 + 6',lM3,12)?^2 

n + 9,2 {{U3,i),2ni - (n3,i),ira2) + 9,i (-(U3,2),2?^l + (ti3,2),l'^2) 
n + 6^,2 ((m3,i),2^2 + (m3,i),1^i) + ^,1 (-(W3,2),2t2 " (M3,2),ltl) 

n + 6',2(%,i),t - 9^i{u3^2),t, 



9,2lU3,2 — 9,i2U3,2 
-9,2lU3,l + ^',12^34 

6',2l'i^3,2 — ^,12^3,2 
-9,2lU3,i + 6',i2M3,l 

9,2lU3,2 — 9,i2U3,2 
-9,2lU3,l + 9,uU3,l 

which is determined by ^3 and n3 „, where we write n = (77,1,71.2)"^. Likewise, 
we can compute that 

div(M(V^ ® Vn3)^) • n 

6',2M3,2),l - (6',2M3,l),2 



',1^3,2), 1 + (6',lM3,l, 



n 



(^^TnU32 12^3^1) ' " (^.2M3,21 - 9,2U3,i2)ni + {-9 ,iU3,21 + 6',lM3,12)n2 



^,21'"3,2 — ^,22^3,1 
-^,11^3,2 + 6',12M3,1 



n. 



which is also determined by U3 and U3,n- Having these computations, we 
conclude that the boundary data {div(LV^?/')-n+(LV^?/')n-t)_t, (LV^Vju-n} 
entirely depends on u' and (^3, 1x3 „) on dVL. 



For the purpose of studying the size estimate problem, we want to show 
that (detV^'?/'), and {hV'^ip) are null-Lagrangians. Throughout the 

paper, (/) denotes the average of the field /, i.e.. 

Using the integration by parts, we immediately obtain that 

l"l Jn l"l Jan l"l Jan 

Hence (V'^ip) is a nuU-Lagrangian. Multiplying the first equation of (2.7) by 
x^/2, X1X2/2, X2/2, respectively, and using integration by parts, it is easy to 
see that (LV^V) is completely determined by the boundary data {^3, u^^n, {h'V'^tp)n- 
n, div(LV2^/;) ■ n + ((LV^V'jn ■ t),t}. 

On the other hand, simple integration by parts implies 



detV = / V',lV',22'^l - V',lV',12'^2 = / V',2V',11'^2 - V',2^,12'^1- (2-13) 

Jan Jan 
Recall that t = — R_Ln. Then we obtain from (2.13) that 

(detVV) = /" (VV)R±n-R^V^ = --i- /" (V2^)t ■ R^V^. 
Jan Jan 

(2.14) 

In view of the expression of Vip on dQ, we have that 

R^V^ = V^,nR±n + ^,tR±t = -i/j^nt + ^,tn. (2.15) 
Substituting (2.15) into (2.14) gives 

(detV'V') = ^ / tn{^'4^)t • t - ^ / ^,t(VV)t ■ n. (2.16) 

I I •/ I I 

Easy computations show that 

(V'^)t ■ t = ^^tt and (V2V^)t ■ n = = (t/^,^),^ 
and (2.16) is equivalent to 

(detVV) = ^ / tni^,tt - ^ / ^,*(^-),* 

I I •/ I I t-* 

l"l Jan l"l 
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In other words, (detV^^) is completely determined by the Dirichlet data 

Our next observation is that 
I 

-divdiv {R{Ve ® Vug) + R(Ve ® Vu^f) ^ + div(MV^VV^)M3 

is a nuU-Lagrangian as well. Indeed, using the integration by parts, one can 
easily check that 

^divdiv {R{Ve ® Vus) + ^(V^ ® Vu-^f) ^ + div(MVVV^)n3 

/ I ( -^div (R(V^ ® VM3) + R(V^ ® Vms)^) ■ n 
l"l Jdn IV ^ 

^ (M(V^ ® Vng) + R{V9 ® Vug)^) n ■ j ^/^ 

^ (M(V^ ® Vug) + R{V9 ® Vug)^) n • tp^n - RV^'(pV9 ■ nug 
= -((V^®Vng) -MV^^A-RVVve- Vug) + fi 

= 5, (2.17) 
where B is the boundary integral term given above. Thus, B is determined 

by V,^,n,%,%,n- 



3 Bounds on the volume fraction 

Now we are ready to derive bounds on the volume fraction of 2-phase shallow 
shell in (1.1) by one set of Cauchy data: 

{(u ; ug, us^n), (s'^ii; (div m - s^V9)n + (mn ■ t),j, mn ■ n)}. (3.1) 

From the above result, we can see that this problem is reduced to the same 
problem for (2.7) with corresponding Cauchy data, which is completely deter- 
mined by (3.1). Moreover, the boundary values of {(ug,Ug^„); (div(MV^M3) • 
n+((MV2ng)n-t),t,MV2ugn-n)} and (div(LV^'0) " n+ ((LV^V')!!- 

t) LV^^/ti- n)} are also determined by (3.1). Therefore, using (2.7), (2.17), 
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and integration by parts, we obtain that (LV^'?/' ■ V'^ip) + (MV^Ms ■ V^Ms) is 
a nuU-Lagrangian. Denote (LV'^ip) = bo, (MV^Ms) = bo, 

A = {(LV^w) = bo] w = ip, w^n = 4',n on dfl}, 

and 

B = {(MV\) =60; V = Us, V^n = Us,n Oil dQ} . 

It is obvious that 

eo := (LV^V^ ■ V^^A) + (MV^s ■ V^s) 

> min(LVV ■ V^^/') + min(MV^M3 ■ V^g). (3-2) 

It should be noted that min(LV^V ' V^V) (resp. min(MV^n3 ■ V^u^)) 

does not necessarily occur at the solution of the corresponding thin plate 
equation divdiv(LV^w) = with Dirichlet data w = ip and w^n. = '^,n (resp. 
divdiv(MV^f ) = with Dirichlet data v = and v^n = u^^n) since we do 
not know whether (LV^w) = 60 (resp. (MV^f ) = 60) holds. To be precise, 
the boundary value problem 

Jdivdiv(LV2u;) = in 

]w = 'iIj, W^n = '4^,n OU dVt 

has a unique solution and thus the corresponding Neumann data {div(LV^ty)- 
n + ((LV^tf )n ■ t).t, (LV^u;)!! ■ n} are completely determined by ip, ip^n- Re- 
call that (LV^tf ) depends on the Neumann data. It is not known whether 
(LV^tf) = (LV^'?/')(= ^o)- However, we have shown that for ip G A (resp. 
M3 e B) we can determine (V^ip) := ao (resp. (V^u^) '■= Oo) and (detV'^ip) := 
Co (resp. (detV^M3) := cq). Therefore, it follows from (3.2) that 

eo - 2Cco - 2Cco > min(LVV ■ V^^^) - 2Cco + min(MV^M3 ■ V\s) - 2C,co 

> min {A ■ LA) - 2Cco + min {B_ ■ MB) - 2C,co, 



(LA)=6o (MB>=6o 
(detA)=co (dctB)=co 



(3.3) 



where C, and C are two parameters which will be specified later. It is clear 
that the two minimization problems on the right hand side of (3.3) can be 
treated separately. 
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We first consider 



min {A ■ LA) - 2Cco. (3.4) 

(A)=ao 
{LA>=6o 
(detA)=co 



Introducing the basis 
V -- 

then any symmetric matrix 



72 VO V'y2V0 -V'V2V1 



A 



Oil 

0/12 0,22 



is represented by a three-dimensional vector v = :^(oii +022, oii — 022, 2012)' 
and the determinant of A is written as 

detA = ■ Tv 

with 

T = 

Similarly, in the basis P, the 4th tensor L is expressed by the 3x3 matrix 





with 

A + 2/i 1 

a = — ; — and p = — . 

4/i(2/i + 3A) 4/i 

Assume the strong convexity condition is satisfied, i.e., /i > and 2/i+3A > 0. 
We can see that both a and /3 are positive. We denote ai, Pi and 02, P2 the 
values in phase 1 and phase 2, respectively. Define a* = min{ai,a2} and 
= min{/3i, /32}. Introduce the translated tensor 

S = L-CT 
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for ( e (— «*), then 5* is also positive. Denote go = bo — (Tclo- Then the 
minimization problem in (3.4) is reduced to 



min {A ■ LA) - 2Cco > min {A ■ SA). 

(A)=ao {A}=ao 
(LA>=6o (SA)=9o 
{detA)=co 



(3.5) 



The minimization problem on the right hand side of (3.5) is exactly the 
one studied in Milton and Nguyen [7] (see (3.8) there). Using their result, if 
A is the minimizer, then 

{A-SA) = go-ao + ^^^^p^[9o-{S)ao]-[if2Si + hS2)go-S^S2ao]. (3.6) 

J1/2 

Applying the same arguments to 



min (B-MB)-2Cco> mm (B ■ SB), 
{B)=do —' ^ - {B)=ao —' 

(SB>=90 



(3.7) 



(detB)=co 

we obtain that if B is the minimizer then 

=^o-ao + ^^^^^|^[^o-(^)So]-[(/2^i + /i^2)^o-^i^2So], (3.^ 



where 



/1/2 

{ S = M -CT, 

a 



C G (-/3=^,a*), = min{ai,a2}, /?* = min{/3i, /^a}, 
. ^0 = ^0 - (Tao. 
Combining (3.3), (3.5), (3.6), (3.7), and (3.8) gives 

Co — 2(co — bo ■ ao + 2Cdetao — 2(co — bo ■ ao + 2Cdetao 



> 



+ 



/1/2 

(^1-^2)-^ 
/1/2 



[90 - {S)aQ] ■ [(/a^i + /iS'2)fi'o - SiS2ao] 
[go - {S)ao] ■ [{f2Si + fiS2)go - 5'i5'2ao]. 
12 



(3.9) 



Expanding the first term on the right hand side of (3.9) reveals 



/1/2 
(5i - 



[go - {S)ao] ■ [{f2Si + fiS2)go - SiS2ao] 
[{bo - (Tao) - {L)ao + (Tao] 



{Si 






/1/2 


(Si 


- ^2)-^ 




/1/2 


{Si 


- s^y' 


/1/2 



/1/2 

■[(/2L1 + /1L2 - CT)(6o - (Tao) - (L1L2 - CTL2 - CLiT + eT^)ao\ 

= ^'^'7f^ [ho - {L)ao] ■ [{f2L, + fMbo - (Tbo - C{f2Li + f^Tao 
J1J2 

~{LiL2-CTL2-CLiT)ao]. (3.10) 
Therefore, this term is hnear in ( with leading coefficient 

-[60 - {L)ao] ■ [(TL2 + L,T)ao - Th^ - (/2L1 + /iL2)Tao] 
-[60 - {L)ao\ ■ [fiLiTao + /2^27'ao - Tbo] 

-[bo-{L)ao]-[{L)Tao-Tbo]. (3.11) 

Note that L2 and T commute. Likewise, we can expand the second term 
on the right hand side of (3.9) and obtain that it is linear in ^ with leading 
coefficient 

i^l^^^[bo - {M)~ao] ■ [{M)T~ao - Tbo]. (3.12) 
J1J2 

Putting (3.9)-(3.12) together, we thus prove that 
Theorem 3.1 The following bounds hold: 
eo-bo- ao-bo- ao 

> C ( ^'^'Tf^ ' [^0 - {L)ao] ■ [{L)Tao - Tbo] + 2co - 2detao 
L /1/2 

+ ' [bo - {L)ao] ■ [{f2L, + hL2)bo - L.L.ao] ^ ^ 

J1J2 (3.13) 

+ ^^ ^^'~fjf ' [^0 - (^)«o] ■ l(M)Tao - Tbo] + 25o - 2detao 
/1/2 
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for all C G (-/?*, a*) and ( e (-/?*, «*). 

We remark that (3.13) is also valid at the extremal points (C? C) = 
a*), or (a*, a*). Therefore, by substituting the 

extreme points into (3.13), we obtain the tightest bounds on the volume frac- 
tion fi (or /2) given by these inequalities. Moreover, when the shell is flat 
{6 = 0) and the plate measurement is either under pure stretching [us = 0) 
or pure bending [ui = U2 = 0) then the bounds derived here reduce to those 
derived in [7]. Note that the pure bending case is equivalent to the pure 
stretching observed in [4]. 
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